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In this paper, we derive equations of motion for the normal- order, the symmetric-order and 
the antinormal-order quantum characteristic functions, applicable for general Hamiltonian systems. 
We do this by utilizing the 'characteristic form' of both quantum states and Hamiltonians. The 
equations of motion we derive here are rather simple in form and in essence, and as such have a 
number of attractive features. As we shall see, our approach enables the descriptions of quantum 
and classical time evolutions in one unified language. It allows for a direct comparison between 
quantum and classical dynamics, providing insight into the relations between quantum and classical 
behavior, while also revealing a smooth transition between quantum and classical time evolutions. 
In particular, the h —> limit of the quantum equations of motion instantly recovers their classical 
counterpart. We also argue that the derived equations may prove to be very useful in numerical 
simulations. 



I. INTRODUCTION 

Quantum characteristic functions play an important 
role in quantum optics, serving as alternative descriptions 
of quantum states - descriptions which are advantageous 
over other equivalent formalisms in certain situations, for 
instance, in the evaluation of expectation values of oper- 
ator moments jlj. Quantum characteristic functions arc 
also very useful in describing the dynamics of open sys- 
tems such as the damped harmonic oscillator, the quan- 
tum Brownian particle, two-level atoms and lasers (for 
reviews, see for example @, S, Q). For these kinds of 
systems, one usually reformulates the appropriate mas- 
ter equation (which involves Hilbert space operators) in 
terms of characteristic functions, to obtain a differential 
Fokkcr-Planck equation, whose methods of solution have 
been very well studied [l[. Recently, quantum charac- 
teristic functions have been shown to be appealing from 
an experimental point of view as well, providing valuable 
tools for highly sensitive measurements [|[. 

While the use of characteristic functions for the de- 
scription of quantum states has a long history, to the best 
of our knowledge, a general equation of motion (EOM) to 
describe their dynamics has not yet been formulated 
This is with the exception of Q , where a general EOM for 
somewhat modified quantum characteristic functions has 
been constructed, using a group-theoretical approach. 

In this paper, we suggest a novel derivation of 
the EOMs for the normal-order, symmetric-order and 
antinormal-order quantum characteristic functions which 
makes use of the characteristic form of both quantum 
states and Hamiltonians. As we shall see, this approach 
provides several attractive features: Firstly, the deriva- 
tion of the equations is both very general and straightfor- 
ward. It leads to EOMs which are rather simple in form; 
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they involve neither Hilbcrt-spacc operators nor infinite 
sums (as opposed to, e.g., the EOM for the Wigner func- 
tion), making them very appealing both for analytical 
calculations as well as for practical purposes such as nu- 
merical simulations. 

Moreover, as we shall see, in their final form the EOMs 
turn out to be very similar to their classical counterpart. 
This property avails insight into the relation between 
quantum and classical time evolutions. Furthermore, the 
formulation of quantum and classical evolutions in a uni- 
fied language enables the evaluation of the classical limit 
in a straightforward manner. These issues are discussed 
in greater detail later on. 

The paper is organized as follows. In Sec. [TT] we re- 
call the definition of the quantum characteristic func- 
tions, and review some of their fundamental properties. 
In Sec. IIII1 we carry out the derivation of the EOMs, and 
in Sec. IIVI the relation of the equations to their classi- 
cal analogue is discussed. We shall conclude with a few 
comments. 



II. THE QUANTUM CHARACTERISTIC 
FUNCTIONS 

Three kinds of quantum characteristic functions have 
been widely used in the literature Q. These are the 
normal-order, symmetric-order, and antinormal-order 
functions, defined, respectively, by the following map- 
pings of density matrices p to complex-valued functions 

B 

C (n) (A,M) = Tr[pe* fa V £a ], (2.1a) 

C (s) (A lA i) = Tr[ / 5e l « a+i « at ], (2.1b) 

C (a) (A,M) - Tr[pe^V« at ]. (2.1c) 



2 



Here, £ is a complex variable defined by the real coordi- 
nates (A, p) as 



^(A-iwM) 



(2.2) 



and the operators a and are annihilation and creation 
operators related to the usual position and momentum 
operators by 



P 



h 
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(d f + a) , 
(a* - a) . 



(2.3) 



The relations between the three characteristic functions 
are given by 

C^(A,/i) =e* lc|a C?W(A,A*) = e®' (\, n) , (2.4) 
and they satisfy accordingly the following bounds 



|CW(A,m)| < e^l\ 
l^ (s) (A,^)| < 1, 
\C {a) (X,p)\ < e"5l«l 2 



(2.5a) 
(2.5b) 
(2.5c) 



It will be useful for our purpose to also recall the in- 
verse of (|2.1a[) . which maps the normal-order character- 
istic functions to density matrices 



1 
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dxdp I dAd M c- l(A * +w) C (n) (A,/i)|a)H. 

(2.6) 



Here, |a) is a coherent state, where the complex variable 
a is defined by the canonical coordinates (x,p) via 



'2hw 



(ujx + ip) . 



(2.7) 



We also note that the normal-order, the symmetric- 
order, and the antinormal-order characteristic functions 
(|2.ip may also serve as the definitions of the Glauber- 
Sudarshan P-representation , the Wigner function [l(| , 
and the Husimi Q quasi- probability distribution [ll| , re- 
spectively, by means of a Fourier transform. 



III. EQUATIONS OF MOTION FOR THE 
QUANTUM CHARACTERISTIC FUNCTIONS 



The time evolution of an arbitrary quantum state p 
whose dynamics is governed by some Hamiltonian H is 
given by the von Neumann equation 



8t r 



pH -Hp 



(3.1) 



In what follows, we rewrite this equation in terms of the 
normal-order characteristic function C^(X 7 p). We shall 
first consider one-dimensional systems; the generalization 
to multiple dimensions as well as the formulation of the 
equation in terms of the symmetric and antinormal char- 
acteristic functions will be discussed later. 

Using Eq. (|2.6p , the left-hand-side of the von Neumann 
equation is simply rewritten as 



d . 

1 

4^2 



(3.2) 



dXdp c -^+^ ^^^ la) (c 



As for the right-hand-side, we make use of the fact that 
quantum observables, specifically Hamiltonians, may also 
be formally represented in characteristic form [3[ . In the 
normal-order case, this representation reads 



H = — [ dxdp [ dXdpe- l{Xx+tip) H^(X 7 fi)\a)(a\ 
1-Kh I ./ 



(3.3) 



with the inverse transformation given by 



2nH {n \X,p) = hTr[He l ^e 1 ^} 



(3.4) 



In terms of C^ n \X, p) and H^ n \X, /z), the product pH in 
the von Neumann equation thus becomes 



(>H = -^j dAid/iidA 2 d/i 2 / dsidpid^dfcC^Ax.pi)^ ^^ 

(3.5) 

I 

Next, we identify the characteristic representation of pH, coherent state diagonal form according to the identity 
by expressing the operator |ai)(a2| in Eq. (|3.5p . in a 
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I "I) ("a | = J d 2 aTr[\ai)(a 2 \A(a,a)]\a)(a\ , (3.6) 
where 

A(a,a) = -jj /" d 2 a'e- Q ' (a ~ at) e a ' (Q - a) . (3.7) 

Integrating over the internal variables, Eq. (|3.5p simpli- 
fies to 



pH = dxdp / dXdpC^(X,p)e- z< - Xx+ ^\a)(a\, 



(3.8) 



with Cp^(A, /i) - the normal-order characteristic function 
of the product /3-ff - given by 

C$(A,/i) = /dA'd/i'C^CA'^O^CA-A',/.-//') 



x e 2 



^j(A / +?cj/j.')(A — A'— iuj(p— fj,')) 



(3.9) 



Similarly, we find that the normal-order characteristic 
function of the product Hp is given by 

C { hI(X,h) = fdX'dp'C^(X',p')H^(X-X',p-p') 



x e 2 < 



-(A —iujfi )(A — A )) 



(3.10) 



Substituting the expressions for pH and Hp into the von 
Neumann equation and equating the integrands on each 
side, we arrive at the final form 



C (n) (A,p) = / dX'dp'K^{X,p,X',p') (3.11) 



d 

at" 

x ffW(A- A', A i-M')C (n) (A',Ai') ; 
with the 'normal-order kernel' iV'^^A, p, A', p') identified 



as 



K^(X,p,X',p') 



(3.12) 



2 e &[V(A-V)WM'(M-M')] sin J(Ap' - p\>) 



Eq. (|3.1ip is the EOM for one-dimensional normal-order 
characteristic functions. 



Of course, in cases where the dynamics of the system 
under consideration is governed by a Hamiltonian con- 
sisting of a linear combination of products of annihila- 
tion and creation operators, Eq. (|3.11[) can be further 
simplified. This is due to the fact the characteristic rep- 
resentation of Hamiltonians of the above form are linear 
combinations of delta functions and derivatives thereof 
Q, so the right-hand-side of p. lip may be immediately 
integrated to produce a linear partial differential equation 
for C^ n > (A, p), to be further solved by the usual methods. 

The EOMs for the symmetric-order C^(X,p) and the 
antinormal C^' (A, p) are readily obtained by making 
use of the relations (|2.4p both for C^(X,p) and for 
H^ (A, p). The resultant equations have the same form 
as Eq. (|3.1ip . but with different kernels. The antinormal- 
order kernel is 



K^(\,p,\',p>) 

_ 2 



(3.13) 



2 e -&[V(A-Y)+«V(M-M')] sin h X p' - p\>) . 



whereas the symmetric-order kernel turns out to be par- 
ticularly simple and reads 



A-M (A, p, A', p') = \ sin \(Xp! - pX) . (3.14) 



The simple form of the symmetric kernel enables fur- 
ther simplification of the 'symmetric-order' EOM (Eq. 
(|3.1ip with («)->(s)). It allows one to Fourier transform 
the EOM and to obtain an EOM for the Wigner function 
W(x,p) - the Fourier transform of C^( A, p). Expanding 
the kernel K^> in a Taylor series in h, the right-hand-side 
of the symmetric EOM becomes 



Em=o b ™ / dX'dp 



m'(A — A') — X'(p — p) 



x H^(X- X',pL-p!)C^{X',p'), 



2m+l 



(3.15) 



with b m = i-iy 



my 



Further expanding the 



(2m + 1)!' 

bracketed term above in a binomial series, and Fourier 
transforming the resultant convolutions term by term, 
the EOM becomes 



oo 2m+l 

E 

m=0 k=Q 



E M-i)* 



2m + 1 



:2(2m+l) 



(gk Q2m+l — k 
g^k dp 2m+l-k W ( X >P) 

I 



2m+l-fe 



0* 



Q x 2m+l-k Q p l 



■H lB \x,p)) , (3.16) 



where H^(x,p) is the Fourier transform of H^(X,p). Eq. ()3.16p is thus a generalized EOM for the Wigner 
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function under general Hamiltonian dynamics. 

This EOM may be written in a compact form with the 
help of the following notation: 

dp dx \9p / \dx^ ) ' (3.17 

where A and B are arbitrary functions of the canonical 
variables. Via this definition, Eq. (|3.16| simplifies to 



whereas the symmetric-order kernel is 



d_ 

Of 



W(x,p) 



(3.18) 



= E^(^fc-fe^) W(x, P )H^,P) 

In the special case where H^ s '(x,p) = p 2 /2 + V(x), Eq. 
(|3.18j) reduces to the well-known equation for the time 
evolution of the Wigner function [13| 



dVjx) d d_ 
dx dp ^ dx 



(3.19) 



^ d 2m+1 V(x) d 2m+1 \ 



m=0 



The generalization of the EOMs derived above 
to N dimensions is straightforward. Now, points 
in phase space are represented by vector pairs 
(x, p) = (xi,- ■ ■ ,xn,Pi,- ■ ■ ,Pn), an d their Fourier 
counterparts by (A, fj) = (Ai, • • • , \n, Mi, • • ■ ,Hn)- In an 
exact analogy with the one-dimensional case, it is an easy 
matter to check that the EOM for the iV-dimensional 
characteristic function C^ n > (A, fi) is 



at 



J dA'd//A (n) (A,/x,A',//) (3.20) 



with 



K( n) (\,fi, \',n') 



(3.21) 



= 2 c& [A'(A-A' )+ - 2 M'(M-M')] sin 5( A/ x' - fiX') , 

while the EOMs for the symmetric-order C^- s \X,fi) and 
the antinormal C' a ^(A, fi) are obtained by making use of 
the multi-dimensional form of relations 

V 2 W/x 2 ) 

i 2 i 2 ■ , 2 \ 



C (n) (A,/x) = c£( A +^> 2 )C (s) (A,At) 

= ci( A2+w > 2 )C* (a) (A,/x), (3.22) 



with similar relations for H^(X, fi). The antinormal- 
order kernel then becomes 



A^(A,/x,A',//) 



(3.23) 



| e -£[A'(A-A') + ->'(/x-/V)] sin |(A/i' - nX') 



2 h 



(A, /x, A', n') = - sin - (A// - /A') . (3.24) 



IV. RELATION TO CLASSICAL DYNAMICS 

Next, we show that the von Neumann equation, in its 
form l|3.20p . enables a direct comparison between quan- 
tum and classical time evolutions. To see this, let us write 
down the corresponding classical EOM, namely the Li- 
ouville equation [l4[ , in terms of classical characteristic 
functions. As with the quantum equation, we start off 
with a general one-dimensional system and generalize to 
multiple dimensions later. 

The time evolution of a statistical distribution P(x,p) 
evolving in time under a Hamiltonian H(x,p) is given by 
the Liouville equation 



dP _ dP dH dP dH 
dt dx dp dp dx 



(4.1) 



In terms of the characteristic functions of P(x,p) and 
i/(x,p), defined by the Fourier transforms 



P(x,p) = J dAd/iC (c) (A,^)e- l(Aa;+w) , (4.2a) 
H(x,p) = y"dAd M H (c) (A,^)e- l(Aa:+w) , (4.2b) 
the Liouville equation becomes 

dCiC) ^ri = J dA'd/ZA^A^AV) (4.3) 
x ff^(A- A',^-^)C (c) (A',m') ; 
where the 'classical kernel' K^(X,/i, A',//) is given by 
A^(A,/i,A',/i') = X/jf-iiX'. (4.4) 

This can be readily verified by substituting (|4.2p into 
(|4.ip . In N dimensions the equation generalizes to 

dCic) (^ri = [dX'dn'KV(X,n,X',n') (4.5) 
x ffW(A- A',/x-//)C (c) (A',/x'), 



with 



A (c) (A,/x, A',//) =X fi'- n- X' 



(4.6) 



Having the quantum EOM (|3.20|) and the classical 
EOM (|4.5p on an equal footing, the difference between 
quantum and classical time evolutions becomes transpar- 
ent. This difference is encapsulated in the dissimilarity 
between the quantum and classical kernels; this is with 
the understanding that in both cases the Hamiltonian 
characteristic function is the same. 
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Moreover, comparing the symmetric-order kernel with 
the classical kernel, we obtain the simple relation 

A-M = % S in -K^. (4.7) 

Relation (|4.7|) reveals a smooth transition between quan- 
tum time evolution and classical evolution, one in which 
h naturally plays the role of a 'quantumness' parameter. 
Specifically, in the h — > limit of K^ s > , we obtain the fully 
classical behavior, namely lim^ = K 1 -^. As one 
would expect, the h — > limit of normal-order and the 
antinormal-order kernels also recover the classical kernel 



V. SUMMARY AND FURTHER REMARKS 

In formulating the von Neumann equation as an 
integro-diffcrcntial equation for each of the three quan- 
tum characteristic functions, we have obtained a 'classi- 
cal' description for the time evolution of quantum sys- 
tems, which further enabled the treatment of quantum 
and classical systems on an equal footing. As a bonus 
we have also obtained a generalized EOM for the Wigncr 
function. In the new formulation, the difference between 
the quantum and classical EOMs narrows down to the 
dissimilarity between the respective kernels. Moreover, 
the correspondence between quantum and classical sys- 
tems arises in this formulation very naturally; a straight- 
forward evaluation of the h — > limit of the quantum 
kernels recovers the classical one. This new formulation 



may thus become very useful when one is interested in 
contrasting the behavior of quantum systems with their 
parallels in the classical world, as analogies can be drawn 
between classical phenomena and phenomena generated 
by quantum dynamics. In fact, a similar strategy has re- 
cently proved to be very helpful in the context of connect- 
ing the no-broadcasting theorem with its classical coun- 
terpart [III]. 

We believe that the EOM presented here, may prove to 
be useful for practical reasons as well. The reason for that 
is twofold. Firstly, the quantum characteristic functions 
are scalar- valued and are ensured to be well-behaved (in 
the case of the symmetric and the antinormal functions - 
also bound); when numerical simulations are concerned, 
these properties are of utmost importance, as they allow 
computations to take place in a rather straightforward 
manner and without having to worry about divergences 
and singularities in the evolving systems. This is in con- 
trast with the Glauber-Sudarshan P-representation Q 
for example, which is ill-defined for some perfectly rea- 
sonable quantum states and is thus very problematic in 
this respect. Secondly, the equation derived here is 'clas- 
sical' in nature; it involves neither Hilbcrt-spacc opera- 
tors nor infinite sums. This property makes the equation 
very convenient for simulation purposes, as this type of 
equations is fairly simple to implement computationally. 
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